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Abstract We report on recent improvements to our non-perturbative calculation of the positronium 
spectrum. Our Hamiltonian is a two-body effective interaction which incorporates one-photon exchange 
terms, but neglects fermion self-energy effects. This effective Hamiltonian is diagonalized numerically 
in a harmonic oscillator basis at strong coupling {a = 0.3) to obtain the mass eigenvalues. We find that 
the mass spectrum compares favorably to the Bohr spectrum of non-relativistic quantum mechanics 
evaluated at this unphysical coupling. 

Keywords Light-Front Dynamics • non-perturbative • bound state • positronium 


1 lutroductiou 

Basis Light-Front Quantization (BLFQ) [I| is a non-perturbative tool for solving bound state problems 
in quantum field theory. This Hamiltonian-based approach combines the advantages of light-front 
dynamics with modern developments in ab initio nuclear structure calculations, such as the 

No-Core Shell Model (NCSM) d; 0. In BLFQ, the quantum field theoretical bound state problem 
is formulated as a large, sparse matrix diagonalization problem. State-of-the-art methods developed 
for NCSM calculations can then be used to address hadronic systems [U;i. The diagonalization 
of the light-front Hamiltonian in a Fock-space basis yields the mass eigenstates of the system, along 
with amplitudes for evaluating non-perturbative observables. Here, we extend previous work Q on the 
mass spectrum of positronium in BLFQ. We evaluate the spectrum at a strong coupling of a = 0.3 
to make both the binding energy and hyperfine splittin g ea sier to resolve numerically. Details of the 
calculations and additional results can be found in Ref. llOll. 


2 Basis Light-Front Quantization 

In BLFQ, hadron observables are evaluated by solving the eigenvalue equation P^Pfj. \'I') = \<F ), 

where P'^ is the energy-momentum 4-vector operator and M is the invariant mass. In BLFQ, we 
express the operator in light-cone gauge. The operator P^ then plays a role analogous to the 
Hamiltonian operator in non-relativistic quantum mechanics. As such, it is sometimes referred to as 
the “light-cone Hamiltonian” H^c = P^ ■ (Note that in this convention the “Hamiltonian” has energy 
squared units.) This operator can be derived from any field theoretical Lagrangian via the Legendre 
transform. In BLFQ, Eq. m is expressed in a truncated basis, and the resulting finite-dimensional 
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matrix is diagonalized numerically. One then examines the trends in observables as the basis truncation 
is relaxed to estimate the results in the infinite matrix (or “continuum”) limit. In BLFQ, three separate 
basis truncations are made. 

First, the number of Fock sectors in the basis is truncated. The basis space for the diagonalization 
in principle includes an infinite set of Fock sectors. For example, the positronium wavefunction could be 

expressed schematically as |e'''e“)pjjyg = a |e+e “)+6 |e+e“ 7 )+c |e+e“ 77 ) + d | 7 ) + /|e+e“e+e“) H-. 

Here we limit ourselves to only the |e“''e“) and |e“''e“ 7 ) sectors. We do not yet make any attempt to 
examine the limit of increasing the number of Fock sectors. 

Secondly, we discretize the longitudinal momentum by putting our system in a longitudinal box of 
length L and applying (anti-)periodic boundary conditions (BCs). Specifically, we choose periodic BCs 
for bosons and anti-periodic BCs for fermions. Thus p'^ = where j is an integer for bosons, or a 
half-integer for fermions. For bosons, we exclude the “zero modes”, i.e. j 7 ^ 0. In the many-body basis, 
all basis states are selected to have the same total longitudinal momentum where the 

sum is over the particles in a particular basis state. We then parameterize using a dimensionless 
variable K = such that P'^ = ‘^K. For a given particle i, the longitudinal momentum fraction 

X is defined as Xi = pf jP^ = ji jK. Due to the positivity of longitudinal momenta on the light-front 
31, fixing K also serves as a Fock space cutoff and makes the number of longitudinal modes finite 
1^. It is easy to see that K determines our “resolution” in the longitudinal direction, and thus our 
resolution on parton distribution functions. The longitudinal continuum limit corresponds to the limit 
L,K ^ 00 . 

Finally, in the transverse direction, we employ a 2D Harmonic Oscillator (HO) basis. The basis 
functions are the eigenfunctions of the operator (q^ = and = y^r^) 
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The eigenfunctions !F™‘(qi) are characterized by the quantum numbers and rrii and the energy 
scale b = VP+J7. The basis is made finite by restricting the number of allowed oscillator quanta in 
each many-body basis state according to (2ni -I- \mi \ -I- 1) < fVmax- The transverse continuum limit 
corresponds to iVmax —>■ 00 . In addition, we use an “M-scheme” basis. That is, our many body states 
have well defined values of the total angular momentum projection Mj = (rrii + Si) , where s = 
is the fermion spin, but they do not have a well-defined total angular momentum J. 

In BLFQ, we construct our many-body basis in single-particle coordinates. The rationale for doing 
this is its straightforward generalization to a basis of many particles. In principle, relative (Jacobi) 
coordinates could be used, but this process rapidly becomes intractable as the particle number is 
increased, due to the need for proper (anti-)symmetrization of the basis states. The disadvantage of 
using single-particle coordinates is that the center-of-mass (CM) motion of the system is contained 
in our solutions. The use of the HO basis combined with the iVmax truncation allows for the exact 
factorization of the wavefunction into “intrinsic” and “CM” components, even within a truncated 
basis. 

When the Hamiltonian is expressed in terms of the coordinates q = p/\/x and s = y/xr, exact CM 
factorization is achieved for all eigenstates, even in a basis with arbitrary numbers of sectors, which 
is the reason for the introduction of these coordinates [HI- The spurious CM excited states can be 
removed from the low-lying spectrum by adding a Lagrange multiplier proportional to Hcm to the 
Hamiltonian to get H' = H + \ [Hcm — 2b'^l) , where H = Hlc- In practice, one selects A to be large 
enough that 2Xb^ is well above the excitation spectrum of interest. Demonstrations of the exact CM 
factorization within BLFQ are given in Refs. [Ill:[l^. 


3 Two-Body Effective Interaction 

We truncate the Fock space to include only |e+e“) and |e+e“ 7 ) states. We wish to formulate an effective 
potential acting only in the |e+e“) space that includes the effects generated by the le+e^y) space. The 
total Hamiltonian can be expressed as (/| iLtot |i) = (/| {Hq + Llmst + Hes) |i), where states |i) and 
I/) are states in the |e+e“) space. The basis states |z) and |/) are eigenstates of the free Hamiltonian 
(i.e. Hq \i) = Ci |i)) with eigenvalue where the sum runs over particles (of mass 
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rrij) in the state |j). The instantaneous photon exchange term, the two-body interaction is the 

only term in the light-cone Hamiltonian directly connecting two states within the |e+e“) sector. We 
choose the Bloch form of the effective Hamiltonian. The Bloch Hamiltonian [T^ is given by: 








( 1 ) 


Since we are interested in primarily the effects of repeated photon exchange, we will only include those 
combinations of terms in Eq. © which generate the photon exchange. We neglect the combinations 
which result in the photon being emitted and absorbed by the same fermion. That is, we do not 
incorporate the fermion self-energy, and therefore no fermion mass renormalization is necessary in 
this model. In addition, we work with unit-normalized eigenstates and a fixed value of the coupling 
constant. 

One can show that iJinst is completely cancelled by a corresponding term in iJeff leaving (a = jA t:) 


.Hinst ^ ^ \J XiX2X-^X2 

aia'^a2^2 


f d'^qi d^q2 (Pc^2 (27r)^j(^)(yTTqi -H - y^q'i - \/4q2) 

J (27r)2 (27r)2 (27r)2 (27r)2 - x[)^[{ei - en) + {ef - en)] 




Cti ,012 fOii ,<^2 


(y^qi, yilq2, \/4q2)> 


( 2 ) 


where we use cti to represent the set of discrete quantum numbers (ji, Si,ni,mi) and to represent the 
subset {ji, Si). The explicit expression for the spinor part, Sai^a 2 ,a[,a' 2 {\/^q.i, v^qzj V^^qii V^q^)’ 
as well as details of the evaluation of the highly oscillatory 8D integration and the derivation of (I2|), 
are presented in Ref. la¬ 
in the energy denominator, we introduce a fictitious photon mass /r to regulate the expected 
Coulomb singularity that, while integrable, introduces numerical difficulties. We then need to examine 
the physical limit p, —>■ 0. 

Previous studies investigating the problem of positronium on the light-front with a one-photon 
exchan ge k ernel have noted a small but noticeable dependence on the ultraviolet cutoff of the theory 
[S [ililS El • Both Refs. [aE| state that the origin of this instability can be traced to a particular 
term in the effective interaction (or one-photon exchange kernel) which tends to a non-zero constant 
(in momentum space) at asymptotically large momentum transfer, corresponding to a Dirac delta 
potential in coordinate space. Since the 2D Dirac delta potential well has no bound states of finite 
binding energy [l^, this leads to a divergence. In our calculation, this would result in a ground state 
energy that does not converge with respect to A^max- If this term is simply dropped, the results become 
convergent. This procedure has some justification, as this divergent piece of the effective interaction 
will be cancelled when higher Fock sectors are included [T 5 j |. Below, we consider the continuum limit 
using only this modified interaction, where convergence can be expected. 


4 Numerical Results 

We now show the results of our numerical diagonalization of the LFQED Hamiltonian using the effective 
two-body interaction described above. Our numerical results were obtained using the Hopper Cray XE6 
and Edison Cray XC30 at NERSC. In principle the basis energy scale b is arbitrary, as any complete 
basis can represent the positronium wavefunction. However, the convergence rate with respect to A^max 
does depend strongly on b. To find the optimal value of b for the ground state, we plotted the ground 
state energy as a function of 6 at AT = A^max = 25, effectively using 6 as a variational parameter, and 
found a minimum at 6 = 0.4m/. This value of b is thus the optimal one for the convergence of the 
ground state energy. 

A representative spectrum of the unmodified interaction is shown in the left panel of Fig. [T] These 
results are produced with K = A^max = 19 and b = g. = 0.1m/. The energies shown are only represen¬ 
tative and should not be considered converged or final. The general features of the spectrum shown 
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Fig. 1 Representative spectrum of positronium {a = 0.3) calculated in BLFQ sd, K = A^max = 19 and 
b — ^ — 0.1m/, using the unmodified (left) and modified (right) effective interactions. The exact energies 
shown should not be interpreted as final converged results. Using the unmodified interaction, the approximate 
rotational invariance allows for the clear identification of the l^So, l^Fi, 2'^5'o, 2®5i, 2^Pi, 2^Po, 2^Pi and 2 ^P 2 
states of the positronium system (see text for details). Using the modified interaction (right), the approximate 
rotational invariance is more strongly broken. 


here are common to any calculation with K = = 19 and above. Convergence will be considered 

below. 

The total angular momentum J of the states can be inferred by examining the multiplet structure 
of the spectrum as it appears in the left panel of Fig. [T] The ground state, for example, appears only in 
the Mj = 0 calculation, suggesting that it has J = 0. We also see a triplet of states above the ground 
state with Mj = —1,0,1, suggesting that these states form a J = 1 multiplet. The lack of manifest 
rotational invariance due to the truncation is seen only in the lack of exact degeneracy between the 
states in this multiplet. The difference, however, is quite small, being approximately 1% of the binding 
energy. We therefore identify the low-lying J = 0 and J = 1 multiplets in the BLFQ spectrum as 
the expected 1^5'o and 1^5'i states of positronium. The remaining states can be identified, via similar 
reasoning, to be the 2^50, 2^5'i, 2^Pi, 2^Po, 2^Pi and 2 ^P 2 states. 

In the right panel of Fig.[T] the spectrum of the modified interaction is shown. The only difference is 
that the rotational invariance is more severely broken. Compared to the basic interaction, the Mj = 0 
states are shifted upwards, while the Mj = ±1 states remain essentially unchanged. This results from 
the ad hoc dropping of a divergent term in the interaction, as discussed in Sec. |31 Despite the lack 
of (near) rotational symmetry in the modified interaction, we will continue to use the same state 
identifications used in the unmodified case. 

We now examine the convergence of these states as a function of our regulators iVmax and K. In the 
left panel of Fig. [U we show the dependence of the singlet and triplet ground states on the regulators 
fVmax and K for a fixed value of the fictitious photon mass fi = 0.04m/. The results are plotted as 
a function of 1/A^max for various fixed values of K as indicated in the figure. The energies converge 
rapidly with K, as one can see from the fact that the curves for K = 45, 55,65 are nearly overlapping 
on this scale. We plot the results as a function of I/iVmax so that the continuum limit (A^max —>■ oo) 
corresponds to 1/A^max = 0. To extrapolate to the continuum limit, we fit the curves to simple second 
order polynomials, although the singlet ground state has a near linear dependence on I/A^max- Tbe 
energies are seen to change by only ~ 0.001m/ over a large range of from A^max = 19 to A^max = 50, 
indicating convergence. 

For a representative higher-excited state, we consider the 2^P2 state, found in our calculation as 
the lowest state of the Mj = 2 sector. In the right panel of Fig. [51 we show this state’s dependence 
on our regulators A^max and K, for the same fixed value of /t = 0.04m/. This state also shows rapid 
convergence with respect to K, however the dependence on 1/A^max is much steeper than for the 
ground state, indicating slower convergence for this state. The reason is that the chosen value of the 
basis energy scale parameter b, which was selected to be optimal for the ground state, is sub-optimal 
for this excited state. In fact, a variational calculation suggests that b ^ 0.1m/ would be optimal for 
this state. 
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Fig. 2 Left panel: Dependence of the and states on the regulators A^max and K for /r = 0.04m/. 
Results are plotted as a function of l/A^max for the indicated values of K. Convergence with K is rapid. The 
continuum limit of A^max —>■ oo corresponds to 1/A^max 0. To extrapolate to the continuum limit, the curves 
are fit to second order polynomials. Right panel: Same for 2^P 2 state. 



Fig. 3 The results of the extrapolations in 1/A^max (as in Fig. [^for /i = 0.04m/) at a series of values of fj.. 
The physical limit corresponds to /r —^ 0. The predictions of non-relativistic quantum mechanics are indicated 
as crosses on the vertical axis. 


To compare our results to the expected Bohr spectrum of positronium, we examine the physical 
limit of ^ — 7 > 0. In Fig. [3l we plot the results of the l/A^max extrapolations (as in Fig. [2] for // = 
0.04m/) at a series of values of /i. For comparison, the predicted energies from non-relativistic quantum 
mechanics are indicated as crosses along the ^ = 0 axis. These predictions include the usual first order 
“fine structure” corrections to the Bohr spectrum of positronium. However, we neglect the correction 
associated with the electron and positron fluctuating into a virtual photon, since we do not include 
the I 7 ) sector in our basis. To examine the limit /r —>■ 0, we fit the resulting curves to second order 
polynomials and extrapolate the results to fi = 0. We find excellent agreement with the expected 
results. It is also interesting to compare our results to the Discretized Light-Cone Quantization (DLCQ) 
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results of Ref. [T^, also using a = 0.3. While their method of calculation is different, their interaction 
incorporates all the same physics and assumptions as ours. In the relevant limit, they obtain a ground 
state mass of 1.97376m/ and a state mass of 1.97708m/, nearly identical to our results. Note 
especially the slight overbinding of the ground state in both methods. 


5 Conclusions and Outlook 

In this work, we have extended our previous results for the mass spectrum of positronium. We have 
truncated the Fock space to include only the |e+e“) and je+e^y) sectors. We further restricted the 
basis to the \e^e~) sector only by developing a two-body effective interaction, incorporating the photon 
exchange effects generated by the |e“''e“ 7 ) sector. However, we neglected fermion-self energy effects 
arising from the |e“'"e“ 7 ) sector. After dropping a divergent term in the two-body effective interaction, 
we obtained converged results in the continuum limit for the I^Sq, l^S'i and 2 ^P 2 states of positronium. 
These results agree well with the predictions of non-relativistic quantum mechanics. 

A straightforward extension of this work is to include a confining interaction between the fermions. 
The model should then be applicable to heavy quarkonia systems. A natural choice for the confining 
potential in BLFQ is quadratic confinement. Such a quadratic potential is motivated by the phenomeno¬ 
logical success of the “soft wall” AdS/QCD model [23;[2l|. The effective interaction implemented here 
would then be interpreted as providing leading-order QCD corrections to the semiclassical approxima¬ 
tion provided by the AdS/QCD model. The BLFQ basis is ideally suited to describe the wavefunctions 
of systems subject to such a confining potential, and we expect convergence with Nmax to be much 
more rapid than in the positronium calculation presented here. 
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